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1. INTRODUCTION 

Power system control is a key role in fulfilling electricity needs. In addition, the increased 
complexity of the load is also a concern [1]. Direct current (DC) motors are included in the category of types 
of motors that are most widely used in industrial environments, household appliances to as supporting 
devices for electronic instrument systems [2]. In a control system, there are several types of control actions, 
including proportional, integral, and derivative control actions. Each of these control actions has certain 
advantages. Proportional—integral—derivative (PID) controller is a combination of the three types of 
controllers. If each of the three types of controllers is independent, the results achieved will not be good 
because each has its own strengths and weaknesses [3]. 

Different controllers are used to control the speed of the DC motor. The most widely used 
controllers are conventional PI and PID controllers. PIDs are widely used in conjunction with DC Motors in 
industrial applications. This control system works by processing calculations based on the control variables 
Kp, Ki, and Kd to achieve the conditions according to the expected setpoint. The PID is able to produce a 
good output response from the DC motor rotational speed [4]. However, its implementation of adjusting PID 
parameters is complex. In recent years there have been many methods for tuning PID parameters. In a simple 
application, a trial-error tuning method is used to adjust the PID value. However, this method is difficult to 
obtain optimal values. So, it is difficult to adjust the parameters, and it takes a long time and also the control 
accuracy is not good [5]. In recent years, researchers have used many artificial intelligent methods to 
optimize the parameters of DC motors, such as particle swarm optimization algorithm [6]—[8], Harris Hawks 
optimization [9], [10], genetic algorithm [11]—[13], firefly algorithm [14]—[16], flower pollination algorithm 
[17]-[19] and neural network [20]—[21]. 
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This paper will present DC motor control using PID which is optimized using the an aquila 
optimizer (AO) algorithm. The AO was introduced by Abualigah [22]. Aquila optimizer (AO) duplicates the 
behavior of aquila in nature during the process of capturing prey. The aquila is one of the most popular birds 
of prey in the northern hemisphere [22]. The contribution of this paper is: 

— The application of the latest and promising metaheuristic methods namely An Aquila Optimizer algorithm 
(AO) to set the PID parameters in DC motors 

— The achievement of the offered method is tested by comparing it with the seagull optimization algorithm 
(SOA), marine predators algorithm, giza pyramids construction (GPC), and chimp optimization algorithm 
(ChOA). 

This paper has an arrangement, namely the second part, which is about the concept of DC Motor and the 

aquila optimizer (AO) method. The third part is the results and discussion. The last part is to draw 

conclusions from the research. 


2. MATERIALS AND METHODS 
2.1. DC motor 

DC motor is controlled by armature and field [23]. The stator and rotor are important parts of a DC 
motor. The non-rotating part of the DC motor is namely the stator. Rotor is the rotating piece. DC motor with 
anchor control uses armature current as of the controlling variable [24]. The stator field is generated by 
permanent magnets or current coils. The motor torque equation is as follows (tm). 


tm(s) = (K,Kyly)Ia(s) = Kmla(s) (1) 
Tm(s) = Kmla(S) (2) 


Where Km is The permeability function of the magnetic material [25]. The armature current (I4) and input 
voltage (V,) have a relationship in the armature circuit the equation is as follows, 


Va(s) = (Ra + La: S). la(s) + ep(s) (3) 

ep(s) = Kya(s) (4) 
where R, and La are Armature resistance and Armature inductance. ep is back electromotive force. 

Tm(s) = T,(s) + Tals) (5) 

T(S) = Jsw(s) + Ba(s) (6) 


Where T; is the torque jointed to the load. tg is Fault torque. J is inertia of the DC motor and B is damping 
friction ratio. The block diagram of a DC motor can be seen in Figure 1. 


Ta (S) 


V.(S) 


Figure 1. DC motor block diagram [26] 


2.2. An aquila optimizer 

Similar to all birds, Aquila has a dark brown coloration and behind the neck is more golden brown. 
Aquila has speed and agility. In addition, the Aquila has strong legs and sharp claws. this supports catching a 
variety of prey. Aquila has been recognized as an adult deer attack. Aquila builds large nests in mountains or 
other high positions. Aquila is one of the most intelligent and skilled hunters after humans. Like population- 
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based algorithms, the AO method begins with a population of candidate solutions (X). The method starts 
stochastically with an upper limit (UB) and a lower limit (LB) [22]. Each iteration will determine 
approximately the optimal solution, which is called the best solution. 


aol |e ae ee 
bere Cae cae 


Xij = rand x (UB; — LB;) + LB; i =1,2....,mjf =1,2,....n (8) 


Where m is the total number of candidate solutions (population), and n is the dimension size of the problem. 
rand is a random number, the j th lower bound is LB}. the j th upper bound of the given problem is UB;. AO 
algorithm methods that simulate the behavior of aquila during hunting can be grouped into four steps: 
— Step 1: Increased exploration (X1) 

In step 1, aquila explores from the sky to determine the area of the search space to determine the 
position of prey. Aquila identifies prey areas and selects the best areas for hunting. 


X(t +1)= Xpest(t) (1 — 5) + (Xu(t) = Xpest(t) x rand) (9) 
Xu(t) = - N X,(t),N =1,2,....dim (10) 


Where, the solution of the next iteration of t is X,(t + 1). It is produced in the initial search method (X;,). 


Xpest(t) Is the best-obtained solution until t the iteration, this describes the estimated point of the prey. The 


parameter to supervisor the increased exploration via the number of iterations is (1 — 5), The points mean 


value of the current solutions linked at ż th iteration is Xy(t). rand is a random value. The dimension size of 
the problem is dim. The population size is N. 
— Step 2: Limited exploration (X3) 

In the second step, the prey has been found with a high level of altitude. In this position, Aquila will 
circle in the clouds, get into position, and prepare to attack prey. At this step, aquila has selected the area of 
prey. Mathematically, the second step can be formulated in (11)-(18). 


Xo (t + 1) = Xpest (t) X Levy (D) + Xp(t) + (y — x) x rand (11) 
Levy (D)=sx“ (12) 
vB 
r(1+p)xsine(7#) 
OS Se = pat (13) 
ret xp 
y =r x cos(0) (14) 
x =r x sin(8) (15) 
r=r,+U+D, (16) 
@=-w xD, +6, (17) 
6b, =% (18) 


2 


Where the completion of the iteration t produced by the second step of the method is X3(t + 1). he 
distribution function of levy flights is Levy (D). the dimension space is (D). Xp(t) is a random solution 
value with a range of 1 to N. s is a fixed constant value with a range up to 0.01. u and v are a random value 
between 0 and 1. ø is a fixed constant value with a range up to 1.5. x and y are used to describe the spiral 
shape in the search. r; is selected a value between 1 and 20 which is used to fix the number of search cycles. 
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U is the variable multiplied by 0.00565. D, is an integer from 1 to the maximum value of the search space 
variable (dim). w is a variable that has a fixed small value multiplied by 0.005. 
— Step 3: Increased exploitation (X3) 

In step 3, Aquila will be in a position of exploitation that is approaching the prey and giving a 
preemptive attack. This behavior can be represented mathematically by the (19). 


X3(t + 1) = (Xpest(t) — Xg (t)) x æ — rand + ((UB — LB) x rand + LB x ô (19) 


Where the exploitation adjustment parameters fixed in this paper with small values (0,1) are a and ô. UB 
indicates the upper limit and LB indicates the lower limit of the given problem. 
— Step 4: Limited exploitation (X,) 

In method 4, aquilla gets closer to the prey. Prey will be attacked by aquilla on the ground. Aquilla 
walk on the ground and take prey. Prey is attacked by aquilla at the last location. Behavioral modeling of 
aquilla in step 4 can be modeled mathematically as in (20)-(23). 


X4(t + 1) = Qf X Xpese(t) — (G1 X X(t) x rand) — G2 x Levy(D) + rand x G, (20) 
_ ,2xrand-1 

Ste (21) 

G,=2xrand-1 (22) 

G, =2x(1-<) (23) 


Where the solution of the iteration generated by the fourth search method (X,) is X,(t + 1). The quality 
function used to balance the search strategy is Qr. All kind of aquila movements used to track prey is Gy. G2 
is a lowering worth from 2 to 0. It is showed the flight incline of the Aquila applied to adhere prey from the 
first spot to the last spot. The current solution at the t-th iteration is X(t). A random point with range between 
0 and 1 is rand. The current iteration is t. The maximum number of iterations is T. The allocation function of 
the flight levy is Levy(D). 


2.3. The proposed AO for tuning Dc motor 

To get adaptive control for dc motors, especially at points such as overshoot, rise-time, and settling 
time. PID parameters are searched and determined using the AO method. Figure 2 is an illustration of the AO 
method in determining the PID parameters on a DC motor. The detail of DC Motor parameters can be seen in 
Table 1. 


Output 
Speed 


K; 
Kp ++ Kas DC Motor 


i 
: : B PID Controller 


Reference 
Speed 


Figure 2. Proposed method diagram 


Table 1. DC motor parameters 


Parameter Value 
Back emf constant (K;) 0.01 N-m/Amp 
Armature resistance (R) 2Q 
Armature inductance (L) 0.25 H 


Mechanical inertia (J) 0.02 Kg/m2 
Friction coefficient (B) 0.05 Nm/rmp 
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3. RESULTS AND DISCUSSION 

AO Algorithm was written and simulated using MATLAB/Simulink on a laptop device with an intel 
i5 (2.2 GHz) processor and 8 GB Ram. Table 2 is the detail of AO variable. To view the potency of the AO- 
PID, it is compared with SOA, MPA, ChOA, and GPC. The global optima function is used to set the 
achievement of the AO method. Figure 3 is the convergence curve. The parameters of the SOA, MPA, GPC, 
ChOA, and AO methods were used to obtain parameters from the PID. Details of the PID parameters of each 
algorithm can be seen in Table 3. The PID data is used to control the DC motor. 

The DC motor speed step response with speed reference 1 for SOA-PID, MPA-PID, ChOA-PID, 
GPC-PID, and AO-PID controllers is shown in Figure 4. Details regarding the step response of SOA-PID, 
MPA-PID, ChOA-PID, GPC-PID and AO-PID can be seen in Table 4. The proposed AO-PID has the best 
reaction step because it has the fastest constancy. Integral total time-weighted square of error (ITSE) and 
integral total weighted absolute value error (ITAE) were used to measure the performance of AO-PID. The 
ITSE and ITAE equations are as follows: 


ITAE = f; t.e(t).dt (24) 
ITSE =|" t.e?(t).dt (25) 
Table 2 Parameter Of AO Tabel 3 The result PID value 
Parameter Value Method P I D 
Solution Number 20 PID 2.8908 9.9239 0.1259 
Maximum Iterations 50 SOA 3.1251 10 0 
Lower Bound 0 MPA 3.18366 10 2.477545 
Upper Bound 10 GPC 2.6571 9.8923 0.1297 
Dim 4 ChOA 3.2061 10 0.1699 
AO 3.18025 10 0 
Objective space s 1.2 


a-o OOo 0-0 ai OOo dl a a d 


=-=- ChOA 
=-0--- MPA 
-= O- -SOA 
sssneneuns GPC 
AO a) 


Amplitude 
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A Onde E a 
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Best score obtained so far 


5 10 15 20 25 
Iteration 


Figure 3. The convergence curve of benchmark Figure 4. Step response with speed reference 1 
function 


Table 4. Comparison result with reference speed 1 
Controller Overshoot Rise Time Settling Time ITSE ITAE 


PID 1.0066 1.811 2.771 0.3069 0.7944 
SOA-PID 1.0037 1.774 2.809 0.2939 0.7633 
MPA-PID 1.0269 2.413 3.15 0.3563 1.140 
GPC-PID 1.0092 1.796 2.711 0.3155 0.8091 

ChOA-PID 1.0045 1.841 2.841 0.2958 0.7828 
AO-PID 1.0032 1.777 2.829 0.2924 0.7631 


The comparison of the ITAE and ITSE with the other controllers can be seen in Table 4. The ITAE 
value of the AO-PID has a value of 0.7631. It is the lowest ITAE value. The ITSE value for the AO-PID 
method is 0.2924. Meanwhile, the highest value of ITSE is owned by the MPA-PID method, namely 0.3563. 
To test the robustness of the proposed method, a test was carried out by changing the speed of the DC motor. 
The reference speed is set to the initial value of 0.8 for 5 seconds. Next, the reference speed is increased to 1 
for 10 seconds. Finally, the reference speed is decreased by 0.5. In Table 5, the overshoot and undershoot 
values can be seen in detail. In the first step, the highest overshoot value is MPA-PID, which is 0.8215. On 
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the other hand, AO-PID method has the lowest value. It is 0.8026. In step 2, the highest overshoot value is 
owned by the MPA-PID method of 1.005. In the second step, the overshoot value of the AO-PID and SOA- 
PID methods differs very slightly, which is 0.0001. In the third step, the reference speed is reduced by 0.5. 
The worst undershoot value is by the MPA-PID method. It is 0.4864. In Table 5, the AO has the same ITSE 
value as the ChOA method, which is 0.3482. The lowest ITSE value is owned by the MPA method is 0.3468. 
ITSE highest score is GPC with a value of 0.3486. The AO and SOA methods have the same ITAE value 
with a value of 0.7387. The highest score of ITAE is owned by the MPA method, which is 0.7398. On the 
other hand, the lowest value of ITAE is the GPC method of 0.7384. 


Table 5. Comparison results by changing the reference speed 
Controller Overshoot Step 1 Overshoot Step2 Undershoot Step3 Rise Time Settling Time ITSE ITAE 


PID 0.8053 1.0014 0.4968 0.854 12.71 0.3483 0.7385 
SOA-PID 0.8029 1.0006 0.4984 0.8 12.77 0.3484 0.7387 
MPA-PID 0.8215 1.005 0.4864 1.139 14.99 0.3468 0.7398 
GPC-PID 0.8074 1.002 0.4954 0.867 12.65 0.3486 0.7384 
ChOA-PID 0.8036 1.0008 0.4979 0.849 12.79 0.3482 0.7388 

AO-PID 0.8026 1.0005 0.4985 0.798 12.79 0.3482 0.7387 


4. CONCLUSION 

DC motor control is a very interesting field due to the rapid development of control methods. Weak 
parameter adjustment will result in dc motor performance. In this study, the Aquila Optimizer Algorithm 
method was proposed to optimize the parameters of the PID. In conclusion, the AO method has optimal 
achievement. The proposed method can reduce the overshoot of the PID by an average of 0.023% and can 
improve the undershoot of the PID by 0.5%. The proposed method, namely AO which is applied to optimize 
PID controllers, has the best ability. 
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